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Abstract 

The derivatives with respect to order for the Bessel functions Jv(x) 
and Yv[x), where v > 0 and x 0 (real or complex), are studied. Rep¬ 
resentations are derived in terms of integrals that involve the products 
pairs of Bessel functions, and in turn series expansions are obtained for 
these integrals. From the new integral representations for dJv{x) /dv and 
dYv(x) jdu, asymptotic approximations involving Airy functions are con¬ 
structed for the case v large, which are uniformly valid for 0 < a; < oo. 


1 Introduction 

Due to their importance in mathematics and physics, Bessel functions have 
been extensively studied in the literature. Their main properties can be found 
by perusing various classic textbooks on the subject ([1], [7], [11], [13], [16]), 
and, most recently, [14, Chap. 10]. The purpose of this paper is to study 
the order derivatives of Bessel functions. The literature is relatively sparse on 
the properties of these functions, with none currently existing on their uniform 
asymptotic properties. We remark that the order derivatives are used in the 
study of the monotonicity with respect to order of Bessel functions, which in 
turn has applications in quantum mechanics [9]. 

For some integral representations see [3] and [12] (the latter for the special 
case V = ±^), and for certain series representations see [4], [8, §8.486(1)], [10, 
§9.4.4] and [15]. We also mention that in [5] modified Bessel functions play a 
role in the study of the order derivatives of Legendre functions 

In this paper we derive new integral representations for the derivatives with 
respect to order of Bessel and Hankel functions (§2), new series expansions for 
the integrals appearing in these representations (§3), and uniform asymptotic 
approximations for the order derivatives of (x) and W (x) (§§4 and 5). 
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2 Integral representations for order derivatives 
of Bessel functions 


Our first result is as follows. 


Proposition 1 For v > 0, |arg(z)| < tt, and z ^ 0 

(^) r Ji it)r j'' it) it) 


du 


= (z) f dt + vttJi, (z) f 

Jo t Jj, 


dt. 


( 1 ) 


If z is complex the paths of integration must lie in C\(—oo,0] (except at the 
end point for the first integral), and with the end point at infinity in the second 
integral on the positive real axis. 


Proof. For convenience, here and throughout we denote 

(z) = dJu (z) /dv, % (z) = dY,, (z) /dv. 


Then, from the v derivative of Bessel’s equation, we have 

rfdfjviz) dJjyiz) / 9 9\ i T / X 

z^—+ z—+ (z^ - ly^) (z) = 2vJ,y (z). (2) 

We then apply the method of variation of parameters, to obtain 

Ju [z) = iznYi, (z) /J" t~^J^ (t) dt + vt:J„ (z) /“ t~^Jv (t) W (t) dt x.,, 
+CuJ,,{z) + d,jY,j{z), 


ior certain constants Ci, and d,^. Recalling the behavior of Bessel functions as 
z —>■ 0''"(see for example, [14, Eqs. (10.7.3, (10.7.4) and (10.15.1)]), we observe 
that the LHS is bounded, whereas the RHS is bounded if and only if d,, = 0, 
which we conclude must be the case. 

To determine Ci, we compare both sides as z —>■ oo. To do so, we use the 
well-known approximations [14, §10.17(i)] 


(f) * <"> =11 - (df °) I 


— sin (z — \vTT — jTt) 


- 1 


O 


and 


(f)1 } 


-I-cos (z - i^TT - jTt) 


- 1 


o 


(4) 


(5) 
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Taking the v derivative of (|3]), we have (formally) for the LHS of ([3]) 


Ju 


/ TT \ 1/2 . 

( 1 

1 \ 

— sm 

, Z — —VTT — 

-TT 

V2z/ ' 

K 2 

4 ; 


( 6 ) 


This can be verified, for example, by a straightforward modification of the saddle 
point method found in [13, Chap. 4, §9.1] to the integral representation 




r* —OO —TTZ 


-oo+7ri 


texp{—zsinh{t) + vt}dt ( jargzj < ) . 


On the other hand, from ([4]) and (O, we find the RHS of ([3]) takes the form 


/ TT \i/2 

J Sin 1 z — — i/TT - tTT 


2 f 1 1 1 ^ 

- cos Z — -I/TT - -TT ] + U 

TTZ J \ 2 4 


.3/2 


In arriving at this we have used [16, p. 405] 

Jo t 2i/ 

Comparing (jH]) with ([7|) we deduce that Cy = 0, and the result ([T|) follows. 
Proposition 2 Under the same conditions as Proposition 1.1 

BY, (z) ,,, r yj (t) 


(7) 

( 8 ) 


f (t) f 

= vttJi, (z) / — - dt — vnYi, (z) / 

J z ^ J z 


^iJj}IllJJldt-\iTJAz). (9) 


dv ' Jz t 

Proof. From [14, Eq. (10.22.6)] we have, for /i ^ ±i/, 
f (z) Pv (z) _ z {C^ (z) (z) — C^+i (z) 'Di, (z)} ^ (z) 'Du (z) ^ ^ 


{p? - 1/2) 




( 10 ) 

where C is some constant. We first choose C = J and V — Y, and we shall select 
C so that the RHS is hnite ds, p ^ v. To do so, note from [14, Eq. (10.5.2)] 


Thus the choice 


C = 


TT [p? — J/2) 


J^+I (z) w {z) - (z) Y^+i (z) = 2 / (ttz) . 


^ z{Ji,+i(z)W(z)-J^(z)w+i(z)} ^ ^ 

t .0 — - T, - 5 -h Oo, 


( 11 ) 


p^ — 1/2 


ensures that the RHS of (ITOl) is indeed finite as /r —>• i/, for any bounded Cq. 
Thus we can assert that 


L 


oo J^j. (t) W {t) _ zY^+i (z) [ (z) - Ju (z) 


P + V 


p + V 


p — V 


p — V 


zYj^ (z) / 1 (z) diz+l i^z) \ Jfi (z) W (z) 


( 12 ) 


P + V 


+ Co. 
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where Cq is now a specific constant which must be determined. To do so, we 
let /r —>■ and (IT^ then becomes 



Ji, (t) Yj, (t) ^ 

t 2v 


dJ„+i (z) 
dv 


Y„ (z) 


dJu (z) 
dv 


Y,+i (z) 


Ju {z) Y^ (z) 
2v 


+Co. 


(13) 

Letting z —>■ oo with v fixed, we find with the aid of (|4]) and that Co = 0. 
Moreover, it is seen that 


r j. (t) Y, (t) 

J. t 


sin (2z — vtt) 
27rz^ 


+ 0 



for each fixed v. 

We next interchange J and Y in (|T^ . to obtain 


(z oo), 


(14) 



Jy (t) Yjy (t) ^ 

t 2v 



jy (z) 


dYy (z) 
dv 


Jy+l {z) 


Jy {z)Yy (Z) 
2v 


(15) 

We can obtain a similar expression by replacing J by K in (1131) . noting that the 
value of Co may differ. In fact, from comparing both sides of m as z —>■ oo 
and referring to (|3]) and ([S]) we find this time that Co = 1/ i^v). Consequently, 
we arrive at 


Yy^t) 


dt = 


2v 


dYy+i{z) dYy{z) 

^y (Z) --- ly+l 


dv 


dv 


z) 


YyHz)-l 

2v 


(16) 


Finally, eliminating dYy+i {x) jdv from (fT^ and (fTBl) . and using (|TT]) . yields m 


Using the behavior of Hankel functions at infinity, and in a similar manner 
to the proof of CD, one can show that the following result holds. Details of the 
proof are left to the reader. 

Lemma 3 For v > 0 and z ^ 0 

^ = -ni \vHi^^ (z) /°° (t) (t) dt 

-vHP (z) (t)} dt-Hi^^z) . 

Likewise, for dul^'* (z) jdv, in interchange the superscripts (1) and (2) 
and replace i by —i throughout. 


3 Series expansions for integrals involving prod¬ 
ucts of Bessel functions 

Here we obtain a number of series expansions for the integrals appearing in 
§2. While power series expansions exist for Jy (z) and Yy (z) themselves (see 
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[14, §10.15], these are only useful for small z, or for special values of i'. Also, 
expansions of the integrals themselves is of importance in scattering theory, as 
noted below. 

Our results appear to be new, and are complementary to certain results ap¬ 
pearing in the literature. In particular, from [16, p. 152], we have the following 
expansion. 


Lemma 4 For positive x and v 

rx j2 1 °° 

/ (18) 
do ^ „=o 

where Eq = 1 and e„ = 2 otherwise. 


This result was generalized in [6]; in that paper it is shown how series of 
this type allows one to write the spherical harmonic component of the Coulomb 
kernel |r — r'|~^ as a numerically stable series. 

The integral in (fT51) diverges when — 0. However, from ([5]), we see that 
(fTSl) can be re-expressed in the form 



4 it) 

t 



1 


OO 

^ ^ ^riJv+n 

n—Q 


(x). 


The integral here now converges when iz = 0 for each positive x, and limit¬ 
ing value of the RHS can therefore be found as ^ 0. In particular, using 

L’Hopital’s rule, we find that 



it) 

t 


OO 

dt = ^ [ EnJn i^) Jn (^)- 

n—0 


(19) 


To obtain a similar result to (TT^ for the integrals appearing in ([T]) and (|9|), we 
shall require the following. 


Lemma 5 For each positive fixed x 


as 12 ^ OO. 


n J, jt) Y, jt) 

Jx t 



( 20 ) 


Proof. Using 


OO 

•'‘•w = (f)''E(-i)‘ 


(U)‘ 




/dr {i2 + k + 1)^ 
and, assuming temporarily that 12 is not an integer, 

Ji, (x) cos (zztt) — J-i, (x) 


Ui/ ix) = 


sin (i27r) 


( 21 ) 


( 22 ) 


we find from (fT^ that (|20l) holds. By continuity the condition 12 not an integer 
can now be relaxed. ■ 
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Proposition 6 For positive x, and v ^ 0,-1, —2, —3, 


r ^ 1 f; (,) (,) + 

F t VTT 2v TT (i/+ n) 

n=0 '' \ ! 

(23) 

where xp {v) is the Psi function [14, %5.2(i)]. For v = —1, —2, —3, • • • f2dl) holds 
with V replaced by \iy\on the RFIS. 

Proof. Firstly we assume v —2, —3, • • •. Then using 


1 


1 


(^) “ r,^P-l (^) r,^P+l (0 


we have 

{Ju+n (t) Yiy+n (t)} = 5 {Jiy+n (t) (t) — Ji,+„+i (t) YJy+n (t)} 

~^2 1 (0 Y,y-\-n (t) Ji/-\-n (t) m- 

On summing both sides from n = 0 to n = yields, via telescoping, 

N 

y {J„+n (t) Y„+n (i)}' = (t) P,.-! (t) + \ Jv-l (t) Y,, (t) 

n—0 

— ^Ju+N+l (t) y^+AT (t) — 5 Jjy+AT (t) Yy+Tv+l (t) ■ 

We plan to let —>• oo. Next we use 

Cp-i (t) = (p/t) Cp (t) + C'p (t), 

to obtain from (IMl) 

N 

y {Ju+n (t) W+n (t)}' = (l^/t) Jv (t) W (t) + 5 {Ju {t) Y„ (t)}' 

n—0 

— {{v + N) /t) Jv + N (t) W+AT (t) — i {Jiz+AT (t) Piy + AT {t)}' ■ 
Integrating this from t = x to t = oo, and using (SD and then yields 

N 

2 ^ 1 / (^) (x) (x) V^_|_ 7 V' (x) (x) (x) 

n—0 

= V /y (t) Yv {t) dt- {iy + N) /y t“Oj.+Ar (t) W+w (t) dt. 
But, for fixed a; > 0, we have from (ED) and ([m 


Ju+n (^) hh+n (^) — 


TT{v + n) 2 tt {y + nf 


O 


(24) 


(25) 


as n —>■ oo. Thus, in order to ensure that the series in (l25ll converges as N —>■ oo, 
we rewrite the equation as 


N r 

^Ju (^) Yi, (x) i^Ju+N (^) (^) 1 Ju+n (^) (^) T 

n=0 I 


1 


= :LiM2kS^dt - 3 £ 


TT (:/ + n) 


TT „=o 1 / + n 


1 -{^+N) /y 


(26) 
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Next, it is well known that [14, Chap. 5] 

^ 1 / 1 \ 

^ —— = tp + N + 1) - Ip {i') = \n{N) - tp {v) + O [ — ) (iV -)> oo). 

V + n V -/V / 

n—0 

(27) 

Moreover, from (1201) . 

in (^)+0(l). m 

Using (l?fl) and (1^ . and then letting —>■ oo in (1^ . yields (1^ . 

Finally, for ^ = —1, —2, —3, • • •, we observe that J^, (t) (t) = J-y (t) Y-^, {t) 

([14, §10.4]), and hence (1^ holds with v replaced by j^j on the RHS. ■ 

In the case ly = 0 we take limits for the RHS of (1^51) . and we obtain 


Jo (t) Fo {t) 


- OO 

‘“-It. 


d 


2 ^ {'C>'+n (^) (^)} 

n=0 


(29) 


1^=0 


In deriving this we used L’Hopital’s rule, along with the relation 
^ - V ''(0 =+ 0(0 ( j ^^ o ). 


We note in passing that if we multiply both sides of 
= 0, we obtain the following new result. 


by u, and then set 


Corollary 7 For positive x 


W I (x) Yn (a;) + — 1 

±1 I 


ln(^a;) +7 

TT 


1 


Jo (a:) Yo {x) ■ 


(30) 


Finally, in order to obtain a complementary representation for t ^Y^ (t) dt, 
one can use the identity 

J-^ (t) F-i. (t) = cos (2i^7r) Jj, (t) W (t) + i sin (2i/7r) { (t) - Y^ (t)} . (31) 


(see [14, (10.4.5)]. Thus, on integrating from t = x to t = oo, we obtain 

pea Y^ (t) „oo Jy (t) . o , fr, ^ r°° 

L = —^dt + 2cot(2j.7r)J^ - - - dt 

^ I'o ^ r°° id) 

—2csc(2i/7r) - - - dt, 


(32) 


provided 2i^ is not an integer. Then one can substitute ([5]), (fTSl) . and (1^ (with 
v replaced by —ly for the third integral in (132^ 1 to obtain the desired result. 

It 2iy = p is an integer we can take limits in (1311) , and in particular use 


lim 

2u^p 


cos (2i/7r) Ji, (t) w (t) - J-v (t) F-i. (t) 


J_ A 

27r dv 


sin (2i^7r) 


{Ju (t) Yu (t)} 


v=p/2 


(-1)" g 

27r dv 


{Jv (t) Yi, (<)} 


v=-pI2 
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As a result, we have 


I 


y2 

oo p/2 


(t) 


dt= r 

Jx 


’Jp/2 


it) 


t 


dt 


T:dv 




Ju (t) Yv (t) dt 


(33) 


4 Uniform asymptotic approximations for J^y (ux) 
and U {i^x) 


We now obtain uniform asymptotic approximations for the order derivatives of 
Jj/ (x) and W (x) for large ly, in terms of the Airy functions Ai (x) and Bi (x). 
For brevity we restrict our consideration to real argument and leading order 
approximations, although extensions to complex argument and expansions are 
feasible from our methods. Explicit error bounds are also available from the 
corresponding ones for the Bessel functions themselves, but again for brevity we 
only include order estimates. 

We note in passing the well-known behavior of Airy functions, given by 


Ai (x) 


exp(-|x3/2) 
2y/Ttx^/'^ ’ 


Bi (x) 


exp (|x3/2) 
0Fxi/4 


{x —7 oo), 


(34) 


and 


Ai (x) 


cos (| |x|^/^ - i^r 


,1/4 


Z, Bi(x)-^ 


sm 1 I \xf^ - 4 


u) 


,1/4 


—oo). 
(35) 


Moreover, Ai (x) and Bi (x) have no zeros for 0 < x < oo. 

We shall employ the corresponding uniform approximations for the Bessel 
functions themselves, a brief description of which is given as follows. Let c = 
—0.36604 • • • be the largest real root of the equation Ai (x) = Bi (x), and define 
the weight function of Ai (x) and Bi (x) by 


E [x) = 1 (— oo < X < c) 

E (x) = {Bi (x) /Ai (x)}^^^ (c < x < oo), 


and the modulus function by 


M (x) = {Ai^ (x) -I- Bi^ (x)}^^^ (— oo < x < c), 

M (x) = {2Ai (x) Bi (x)}^^^ (c < x < oo). 


Next, introduce a new variable C by 


= In 


1 -I- (l - X^) 


1/2 ■ 


X 


(l-x^)^^^ (0<X<1), (36) 



- (—— sec ^ (x) (1 < x < oo), (37) 

O 

all functions taking their principal values, with C = oo,0, —oo, corresponding to 
X = 0,1, oo, respectively. Note C —oo as x —>■ oo such that 

x=l{-(:f^ + l^ + o(\C\-^/^), (38) 

and C —>■ oo as X —)► 0+ such that 

X = 2 exp - l) + O |exp | • (39) 

For complex x and C, the transformation is given by (13611 . where the branches 
take their principal values when x G (0,1) and C S (0,oo), and are continuous 
elsewhere. 

From [13, Chap. 11, §10] we then have 


and 


where 


and 


Y,,{iyx) = -^ (-^) {Bi(i/ 2 / 3 c) +ei(i/,C)}, 

(tZ^) {Ai(i/'/3c)+e2(^^,C)}, 


£1 {u, 0 = 0 {u-^) E M , 

£2 {i'. 0=0 {v-^) M ( 1 . 2 / 3 ^) /E , 


uniformly for 0 < x < 00 (i.e. —00 < C < 00 ). 

We shall also utilize the uniform approximation 


(zzx) 




(40) 


(41) 

(42) 

(43) 


(44) 


where (for our purposes) x and ( are complex, with —tt < arg(C) < — ^tt, and 
correspondingly x lying in a certain unbounded subset of the upper half plane 
(see [13, Chap. 11, §10]). A fundamental property of Ai {O'^''l^v^O0^ is that for 
large v it is exponentially small in the sector — tt < arg (C) < — ^tt (and likewise 

{vx) for X in the upper half plane and lying outside a certain domain which 
contains the origin: see Figs. 10.1 and 10.2 of [13, Chap. 11]). 

Now, from o and m we have 


Ji, (vx) = vkJv {vx) Ii {v, vx) + vttYi, {vx) I 2 {v, vx), (45) 


and 


% {vx) 


vttJi, {vx) I3 {v, vx) — zzttIA {vx) Ii {v, vx) 


-ttJ,, {vx) , 


(46) 
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where 


{v, t) = 

r JAs)YAs)^^ 

(47) 

J 

't s 

h t) 

II 

.'ll 

(48) 

JO s 

h {y, t) ■- 

Jt s 

(49) 


The Bessel functions in (HU) and can immediately be replaced by their 
Airy function approximations, given above. Regarding the integrals Ij {v, t) 
{j = 1,2,3) we note, from (l40l) - ((43ll and ((35l) . that the integrands in (ITf)) 
and (j49ll are highly oscillatory in the interval v < s < oo^ and likewise for 
(j48ll if t > V. This is detrimental in numerical evaluations, so our goal is to 
approximate these integrals by readily computable functions and/or integrals 
with monotonic or slowly varying integrands. 

Our results are summarized in Theorem 4.2 below. In this section and the 
next, C (t) appearing in an integral denotes the function given by (l36l) and (1371) 
with X replaced by t, and similarly x (ji) is given by (I5S1) and (1371) with replaced 
by V- 

Before stating our main theorem, we shall require the following result. 


Lemma 8 Let 


Then as v ^ oo 


L{v,t) 



(50) 


L Iv, vx) 


4g-2W3 f Cit) \ Ai^ ^ 



(51) 

uniformly for 1 < x < oo. Here the path of integration lies in the first quadrant 
and is bounded away from the boundary EPB in Fig. 10.1 of [13, Chap. 11, 
§ 10], except for the endpoint when it is located at x = 1. 


Proof. We make a simple change of variable s = vt to obtain from dsni) 


L [v, vx) 



{vt)] 

2t 


2 


dt. 


(52) 


We then use (H^ in the integrand, and deform the path of integration from 
the positive cc-axis to the stated path. The deformation is justifiable since 
the integrand is holomorphic in the first quadrant, and exponentially small at 
X = zoo. ■ 

Note that +oo as x —>■ zoo, and hence the Airy function in the 

integrand of m is non-zero and (in absolute value) rapidly decreasing along 
the path of integration. 
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Theorem 9 For large v, uniformly for 0<x<l (0<C< c»), 


Ii (y, vx) = Im L {i', v) 

2 Ai (1/2/3^ it)) Bi {t)) 

^- 


y/3 


t 


dti 1 + 0 


(53) 




and 

h 


h r (B) {i+o (i; 


jy2/3 

For large v, uniformly for 1 < a: < c» (—oo < C ^ 0) 

Ii {v, vx) = Im L (v, vx ), 

I 2 (v, vx) = - KeL{v, vx) 

2v 


(55) 

(56) 


1 


L 


2/3 Jx \l_-j-2 


c(i) A^{vy^c(.t)) + Bi^vy\{t)) 


dt <1 + O 


(57) 


and 


I 3 {v, vx) = — Re L {v, vx) 

1 .no f C{t)\ Ai2 (^.2/3^ (^2/3^ (i)) 


”*",,2/3 /a 


1-^2 


t 


dt <1 + O 


In these, L {v, vx) is approximated by m- 


(58) 


Remark. The Airy functions in the integrals appearing in (l53l) . (l54l) . and 
the first integral on the RHS of (1551) are all non-oscillatory, since their arguments 
are nonnegative. From (1551) we also note that 

Ai2 (j. 2/3^^ ^ Bi^ (z. 2/3^^ _ |^|l/2^ ^ ^ 

and indeed, as seen (for example) from [14, Eq. (9.11.4)], this function is mono¬ 
tonic. Therefore the second integral on the RHS of (l55|) . as well as integrals 
appearing in (EB and (1581) are slowly varying, and consequently all the inte¬ 
grals in question are numerically stable. 

Proof. Making the change of variable s = vt in (|T7)) we have 


h 



Jy (vt) Yy {vt) 
- dt 


(59) 
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(with similar representations used for I 2 (v,vx) and [v,vx)). We next use 
(t) = Jv {t) + iYi, (t) to give the relation 




for t positive. Then, using this in (1551) and referring to ((55|) . we see that (|56p 
follows. 

For 0 < cc < 1 we express Ii {v, vx) in the form 


{u, ux) = / 

J X 


Jy {vt) W {vt) 


dt - 


Ji, [vt) W {vt) 


dt. 


(60) 


The first integral here has a desired monotonic integrand, and using (HOI) and 
m it can be approximated by 


Jx 


Jy {ut) Yj, {yt) 


t 


dt 


1 f Cit) \ Ai (z/2/3( (t)) Bi (t)) 


. 2/3 


£(t 


-t2 


dt < 1 + O 


For the second integral, we have from (l56l) . 

/■°° J, {l2t) W (i^t) 


dt = Ii (y, v) = Im L {v, v) . 


(61) 

(62) 


Thus (l60l) . ((6T|) and (l62)) establishes (l53|) . 

Similarly to (|6T|) . we use (HD) in (|48l) to obtain (l54l) for the case 0 < a; < 1. 
For 1 < a; < 00 we use Q to rewrite 


I2 {v, ^^) = ^ - 




dt. 


(63) 


Next, in this we use 

Jl it) = ^ Re I (h(i) (t))'| + (t) hP it ), 


to arrive at 


■im.t =R. r 




2t 


2t 


(64) 


Thus (1571) follows from (HOI) . (HT1) . (1501) . (100)) and (100)) . 
Similarly to (1001) , we use (HOI) along with 


it) = it) it) - i Re I (t))'|, 


12 













to arrive at (1551) . 

Finally, for 0 < x < 1, we write 



The approximation (I55|) then follows from using (HOI) in the first of these inte¬ 
grals, and using (1^ (with x = 1) for the second integral. ■ 

As an illustration, let v = 100 and x = 0 . 5 . The exact values that follow 
are obtained using the fdijj routine in MAPLE with sufficient precision. From 
(H51) we obtain the approximation Jioo ( 50 ) ~ — 1.47735 x 10 “^^, compared 
to the exact value Jioo ( 50 ) = — 1 . 47702 ... x 10 “^^. Similarly, from (H51) . we 
obtain the approximation lioo ( 50 ) ~ — 4.31473 x 10 ^®, compared to the value 
Yioo ( 50 ) = - 4.31569 x 10 ^®. 

Likewise, if we choose v = 100 and x = 5, we obtain the approximations 
Tioo (500) 

Ri 0.0150731 and Fioo (500) ~ —0.0470087. In comparison, the exact values are 
Jioo (500) = 0.0150695 • • • and Tioo (500) = -0.0470099 • • •. 


5 Asymptotic approximations of the integrals 
appearing in Theorem 4.2 

The integrals in Theorem 4.2 involve Airy functions, and as such it is desirable 
to have explicit asymptotic approximations for them. These will involve the 
following integrals 

Gi,] (y, 0 = j C'"y^ Vj dc, (65) 

for i = 1,2, j = 1, 2, and /c = 0,1, 2; here yi (t) = Ai (t) and 2/2 (t) = Bi (t). The 
key is that these integrals can be explicitly evaluated, and in particular from [2] 
we have (with i and j not necessarily distinct) 

Glj (y C) = Cy^ yj y- , (66) 

GG (y C) = ley. y' .,7. 

{y. y' + y' y, ^ ’ 

and 

GL(-,C) = -i--^/^C^y'(-^/^C)y'(-^/^C) 

+ i^- 4 / 3 C {y, (z. 2 / 3 ^) (^2/3^) + y. (^2/3^) y^ (^2/3^) } ( 68 ) 

Next let 

ao = 2 oi = -, 02 = —2^/®, (69) 
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these being the first three coefficients in the Maclaurin expansion 

= PO) 

see [13, Chap. 11, §10]. Finally, let a:o = 1 — (correspondingly Co = C (a^o) 
= (2/:/)^^^ + O and xi = 1 + (correspondingly Ci = C (a^i) = 

The desired results are stated as follows. 


Theorem 10 For large v, uniformly for 0<a:<l (t)<C<oo) 
C (t) \ Ai (i/2/3C (t)) Bi (t)) 


ri (at) y 

Jx \i-ty 


t 


-dt 


= 1 ^ ^ ) r S “in{c, Co}) - ^ o)] ( 7 i) 


+ 


1 


2'kv^/^ 


1 + O ( — ) ]• max \ cosh ( — ) — cosh ( — ) > 0 ]■ , 


1 


r( 


C(t) AP(^2/3^(t)) 


1-^2 


dt = — 


1 — 


GlAv,0{l + O 


and 


A f c(t) y/^ Bi^(^2/3^(t)) 
■ ' C ' 


1 — 


GL(^,C)-aoG° (z.,0) 


1 + 0 


(72) 

(73) 


In jyJI ) the nonnegative branch of the inverse hyperbolic cosine is taken. 
For large v, uniformly for 1 < x < oo (—oo < C < 0) 


4e-2W3 f ^ 

^ = —f2jT 


1 — 


G?.i (e"*/3^^,c) 


1 + 0 


(74) 


and 


C (<) \ Ai^ (^.2/3^ A)) + Bi^ (^.2/3^ A)) 


l-t2 


dt 


Ffi 


1 + 0 ( - 
V 


sin 


1 


max {x, xi} 


+ max \ J2 Ofc [Gyi A, C) + G^ 2 C) - Gf_i A, Ci) - G§ 2 A, Ci)]. 0 




(75) 


with the principal inverse sine taken in m- 

Proof. For the integrals (17^ . (1751) . and (1771) that have rapidly changing inte- 
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grands (i.e. exponentially decreasing or increasing), we use integration by parts. 
We consider the latter, with other two results being similarly proven. Applying 
integration by parts to the integral appearing in (1511) we obtain 






Ai2 ^g27ri/3j^2/3j^^ _ 


/ooe-2-i 


/3 


1 — a; {riY 


1 — a;^ 




(76) 

we 


Here the prime represents differentiation with respect to 77 , and from 
have 

Gi.i = CAi^ • (77) 

Now the first of (j34p holds for complex x with |arg (a;)| < tt — i5 (J > 0) and the 
principal value of taken. Then, from that equation, along with (16511 and 
dZZl), we can readily show that G? ^ Y) is non-vanishing for nonnegative 

C, with 

6 - 2 -/ 331 / 3^2 ^ 2 ^ 


G' 


’ 1 


dTT^ 


^ -+ 0 (C). 


as C —t 0 , and 


G' 


'1 (e"/"^,C) 


6 exp {|z^C^/^} 

87rz/4/3^ 


1 + 0 


l/3^3/2 


as 


j^^ 3/2 ^ such that |arg | < tt — S. Hence we can assert that 


K, 


|exp{|K"/ 2 }| 
„ 4/3 1^1 + ^ 2/3 


< 


G°i (e-/3za,c) 


< K2- 


I exp I 

l/"‘/3 |()| + zy2/3 ’ 


(78) 


for some positive constants Ki and K 2 . Next, from the complex variable exten¬ 
sion of (l36l) it is straightforward to show that 


c 


1 — X" 




^ I ^3 / > 


(79) 


as C —t 00 in the sector |arg(—< 27r/3. Therefore from this (with x, ( 
replaced by x{ri), 77 , respectively), and the upper bound in (1751) . we have for 
large u 


Jooe“373i/3 ^ 

= 0(l)/i 


77 


l-x( 77 ) 

.-27vi/3 


G]* 1 77 ) drj 

I exp 


(1/4/3 1^1 _|_ jy2/3) _|_ X)- 


'd\v\, 
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uniformly for —oo < C ^ 0- Now express r] = on the contour, where 

r = |? 7 |, and 9 ir) decreases continuously from tt to 27r/3 as r increases from |C| 
to 00 . From Laplace’s method ([13, Chap. 3, §7]), we arrive at 


/ - 
J ooe 


27rt/3 


l-a;(? 7 ) 


G° 1 r?) dr] 


= 0{l)l 


ICI 


exp cos (f (?"))} 


dr 


= O 


( i / 2/3 1^1 + 1 ) ( 1^1 + 1)3 


(80) 


The result (1711) follows from (ISTl) . (1751) . (1751) . (17^ and (1501) . The proofs of (1701) 
and (1751) are similar. 

Consider next the integrals with slowly varying (non-exponential) integrands. 
The method of integration by parts is not effective in approximating these, and 
instead we resort to a combination of a Maclaurin series expansion for small 
values of the integration variable, coupled with asymptotic approximations of 
Airy functions for the other unbounded values. To this end, on using 


^ _ _1 

dx X \ ^ / 


we express m in the form 

ji <{t) y/^ Ai(^^/3c(t))Bi(^/^/3c(t)) ^^ 

= fo” ( , \ y ) Ai Bi (81) 

\l-x{r]) J 

^ j.o Ai {y^/^C jt)) Bi (yy^gt)) 

If 0 < C < Co ( 2^0 < a; < 1 ) the upper limit in the first integral on the RHS is 
replaced by C and the second integral is null. 

For the first integral on the RHS of (15T1) we find from (1701) . and recalling 
Co = ( 2 / 1 /)^^^ + O [y~‘^/^), that 


\l-x{ri) I 

= |l -b O I fg° {ao -b ai?7 -b 02 ^} Ai {y'^/^'q) Bi {y'^/^q) dq. 


(82) 


Explicit integration and referring to (lOSl) yields the first term on the RHS of 
GB). Clearly this also holds if Co is replaced by C when 0 < C < Co- 
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Next, from (131 


Ai 


Bi {i 


^ 2 / 3 ^ = _ I _ 

V 27rz/i/3^i/2 


1 + 0 


y2(^ 


as ^2/3^ 

—>■ oo. Hence, assuming 0 < a; < xq (Co < C < oo)) we have C ^ = 
O in the second integral on the RHS of (IHTl) . and so it follows that 


r^o f C {t) \ 
Jx 


C{t) y/" Ai (1/2/3^ (i)) Bi (j. 2/3^ (t)) 


2tti^03 


r( 


1 \‘'"l 


1 -y t 


■dt 

dt. 


(83) 


Now for 0 < t < 1 


1 y/2 1 /I 

- -^ I -dt = — cosh“^ f - 

1 -12 y t 


(84) 


On combining ((Ml) and (|M|) we then arrive at the second term on the RHS of 

dZl-dMI) 

The proof of (1751) is similar. For 1 < x < Xi (Ci < C ^ 0) we write this as 


c (t) \ Ai^ it)) + Bi" (j^"/3C it)) 


dt 


^ Ai" K=cw)+Bi"K=c(t))j , 


+ /( 


XI _ ^2 

C I V 


(85) 


11 — X iY 


|Ai" iy^Or}) + Bi" drj, 


with no such splitting being necessary if —oo < C ^ Ci i^i ^ x < oo). Now 
using 


1 




Ai"y2/3y+Bi"(z."/3c)=- ^ 

V V V y (_^)i/2 

we find for the first integral on the RHS of ((85ll that 

C(t) V/"Ai"y/3^(t))+Bi"y/3^(t)) 


roof at) y 

Jxi j 


t 

1 \^/"i 


dt 


Then with 


TTi^Os + Vt 2 -iy t 


1 \^/"l 


-dt. 


-dt = sin ^ [ — ) , 

t2 - ly t Va;i 
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we arrive at the first term on the RHS of GSD. 

In the case when splitting is necessary (i.e. 1 < a; < xi), we follow (l82|) 
and employ the Maclaurin expansion (1701) for the second integral on the RHS 
of (I85|) . and the second term on the RHS of GS) is obtained. ■ 


Table 1 illustrates the results (dD, (EH), (EH and (1751) numerically for the 
case = 50 and various values of x. In this, we denote the relative errors 
m {I = 1,2, 3,4) by 

fi x) = gi {v, x){l + rii {y, x)} , 


where 


and 


9i (v, x)= ak [Gn 2 min {C, Co}) - ^ {v, 0)] 

fe =0 


/2 = 


max < cosh 


— cosh 


-1 


,0 , 


c(t) Ai^juv^cit)) 


i-tV 


dt. 


92 {v, x) = - G^i {v, C), 


53 [v, x) = 


c 


1 -X 2 

1/2 Ai2 


G°2 (z/,C)-aoG°,2(z.,0) 


.. ^ r (at) \ Ai^ a)+ Bi^ a) 
a a X) - - 


dt, 


54 a x) = ^73 sin , - - -^ 

\max|x,xi| 

+ max I X) Ofc [Gyi {v, () + G^ C) - G^.i Ci) - G^ a, Ci)], 0 }. 
U=o 


All calculations were performed with MAPLE, with the integrals // (z^, x) being 
evaluated numerically using Simpson’s method. 


X 

0.1 

\m (50,x)| 
1.6240E;-04 

152 (50,x)| 
3.1202A-03 

0.5 

3.4440E; - 04 

6.9778A- 03 

0.75 

2.1710E;-04 

1.0326A- 02 

0.99 

1 

5 

1.7825E-08 

2.0756A- 02 

_ 

_ 

10 

- 

- 


|773 (50,x)| |54(50,x)| 

3.1466E;-03 

7.2109E;-03 

1.1859E;-02 

1.8467E;-02 

2.6086E-04 

6.2709A-07 

1.1871E;-07 


18 

















Table 1. 
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